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Abstract 



An algorithm for computing 7r(N) is presented.lt is shown that using a sym- 
metry of odd composites we can easily compute 7r(N).This method relies on the 
fact that counting the number of odd composites not exceeding N suffices to 
calculate 7r(N). 
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Introduction 



The problem of finding a compact mathematical formula for 7r(N) is well known 
[1]. Though no compact result has been reported till date, a number of accu- 
rate/fairly accurato algorithms exist [2,3,7,8] among which the prime number 
theorem [4] and the function Li(x) [3] require special mention. Wilson's theorem 
[5] seems to be unique, however, hke all reportings so far it has a low practical 
utility. Recently a result has been reported using Smarandache Functions[6] 
which too however has a low practical utility. 

In the present paper we devclope a computational algorithm based on a cer- 
tain symmetry of odd composites. This is the main essence of the present piece 
of work. 

The Algorithm 

A. General Formulation 

Among all natural numbers we know that 

i) All even numbers except 2 are composite. 

ii) An odd composites can be factorised into odd primes only. 

iii) Unity is neither prime nor composite. 

Therefore we have 

7r(N) = N - ne(c) - no(c) - 1 (1) 

where 

ne(c) and no(c) give respectively the number of even and odd composite 
numbers not exceeding N. Without the loss of generality we assume N>2. 

Now we have 

ne(c) = Floor (f - 1 ) (la) 

To obtain no(c) let us refer to the multiphcation table of odd natural num- 
bers >3 as shown in Table-II. It shows a very small section at the top left hand 
corner of the infinite multipliction table. 

It can be shown that all odd composites occur at least once within Table - 
II whose general term is given by 

tn^r = ( 2n + 1 )2+ 2 ( r - 1 ) ( 2 n + 1 ) = ( 2 n + r )2- ( r - 1 f (2) 

where both n and r are natural. 

We further see that for each 'n' value we have an A.P originating at ( 2n + 1)^ 
with a common difference 2 ( 2 n + 1 ). The terms of the A.P for each 'n' are 
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assumed to constitute a set represented as 
{ t„,, } = { X : X = (2 n + r )2- ( r - 1 )2} 
where r varies from 1 to 

rUax = Floor i^^^i N - ( 2n+l)2) + 1 ) (3a) 

for a given n i.e for a given set { tn,r }■ 

Thus, 

no( c ) = E Floor + 1 ) - (4a) 

where the sum is over n from 1 to 

nlmaa: = Floor (^^) (3b). 

Here Ac is a correction term which we would discuss shortly. 
B. Corrections 



Bl. The Summation 



For a composite (2nr+l),all the terms arising from the summation in (4a) cor- 
responding to n=nf. would be a duplication of some previously occuring term 
for some n = n,, <: n,. .Therefore we shall restrict ourselves to only those values 
of n,- for which (2nr. + 1) is prime. 

Following the discussions as above we can recast the summation part of (4a) as 
no(c) = Ef (Floor ^-g:;+y +l)-A,...(4b) 

where we restrict the sum for a prime (2ni+l). Here the suffix R denotes the 
sum over the restricted sets of the values of n's.The primality of (2nj+l) can be 
ascertained using the symmetry we mentioned. We do not accept those values of 
n/j for which 

nfc <nmax for r < 

^max cis given by (3a) and (3b). 
As win be shown in a latter article [ 9 ] ah composite (2nj + 1 ) can be expressed 
as a difference of two squares following eqn. (2) . 

Alternatively the primality could be determined using Wilson's theorem[5]. 

B2. Evaluation of Ac 

The need of Ac arises from the fact that for different n values the sets { t„^r} 
coincide at certain points. 

Now once the summation in (4a) be restricted we shall also require to restrict 
the evaluation of Ac. 

Restricting ourselves to those values of n^ for which (2nr+l) is prime we can 
evaluate Ac using the set theoretic formula 
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n(AiUA2U--UAfe) = E n(A,)-EEn(A.nA,)+E n(A,nA,nA;c) - (5) 

where n(X) denotes the cardinahty of set X. 

Here summations over all indices can range from 1 to n^ax with the restriction 

i < j < k ... and so on. 

Let us define the following notations: 

For intersoction of two A.P's let 

tni,n2 be the first common term of { tni,r} and { tn2,r } while 
( c.d )ni,n2 be the common difference of the common terms 
For intersection of more than two A.P's the notation is identical. 
It can be shown that 

a ) For two A.P's 

t^nln2= (2ni + l )* ( 2n2 + 1 ) * { 3 + 2 * Ceil [Z]*U[Z]} ... (6a) 
where 

7 _ 1 / 2^2 + 1 q A 

^ - 2i2;iT+T - ) 

and 

(c.d )nun2= 2 * ( 2ni+ 1 ) * ( 2 n2 + 1 ) ... (6b) 

Here H [ Z ] denotes the Heaviside Step Function defined by 



H[z]= 1 , z > 
= , z < 

b ) For more than two A.P's 
ti'ln2...n,= nil ( 2n,+ 1 ) * ( 1 + 2 * Ceil ( Q) * H [ Q ] ) ...(7a) 
where 

and 

( c.d U,n2,...,nk = 2 * ni=i( 2ni + 1 ) ...(7b) 
Using (5), (6a), (6b), (7a) and (7b) we get 
Ac=E;rr(-l)''E'' E''-E''(Floor(B) + l)*H[B+£] ...(8) 
where 

i) The summations run over all restricted sets of n values, one for each set with 
ni < n2 < na < nfe 

ii) e is a positive number, however small. It will be neglected at the last stages. 

iii) We have 
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(c.<i)„j,„2,.,.,„^ 

The Heaviside function restricts the numerator to non-negative values. 

Using (1), (la), (4b), (5), and (8) the result follows. 
Though equation (8) seems to have many steps in practice it would be quite 
Hmited as once the result given by (6a) and (7a) exceeds N, higher n^ terms 
would be in general very very limited in number. 

The Results 

The results as obtained from the above mentioned algorithm are given in Table 

I. It is seen that the value of 7r(N) is exactly reproduced. This estabHshes the 

correctness of the above mentioned algorithm. 

To illustrate the method let us evaluate 7r(100).We have 

from (3b) : nmax=^ 

Now, to calculate 7r(100) we need to know all primes till a/100=10. 
Corresponding to n=4 we get (2n + 1)=9 < 10 which is a composite. 
So we do not consider n=4. 
Proceeding: 

from ( la ): ne(c)=Floor(i§2 - 1) = 49 
from (4b):no(c)= 16+8+4=28 
from ( 8 ):Ac=2+l=3 

Therefore from (1) we have 7r(100) = 100 - 49 - 28 + 3 - 1 = 25 
Conclusion 

The formula for 7r(N) as established seems comphcated at first sight. However 
the symmetry of natural numbers used, the necessity of introduction of Ac 
and the simple set-theoretic prescription for its evaluation makes the present 
method quite tractable. However, like the well known "Sieve algorithm" [8] this 
method too requires the knowledge of all primes till \fN to evaluate 7r(A^) .The 
number of steps is much reduced as we use the symmetry of natural numbers 
throughout, first to determine primality of (2nj+l) if the same is not already 
known and afterwards for the final evaluation of 7r(N) in comparison to the 
Sieve Algorithm. Refering back to the symmetry of odd composites that we had 
mentioned, the symmetry is more clearly evident if we use the binary base of 
representation specially in Table H.This is evident from the fact that the com- 
mon difference of the common differences in Table H equals base squared. It is 
believed that this could be a suitable method to determine 7r(N). 
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TABLE : I 

Results of determination of 7r(N) 



N 


7r(N) 


Estimated A'alue 


10 


4 


4 


100 


25 


25 


1000 


168 


168 



TABLE : II 
Multiplication Table of odd natural numbers 





3 


5 


7 


9 


11 


Common Difference 


3 


9 


15 


21 


27 


33 


6 


5 




25 


35 


45 


55 


10 


7 






49 


63 


77 


14 


Common Difference 












4 =2^ 
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